Extended Fluctuation-Dissipation Theorem for Soft Matter in Stationary Flow 
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For soft matter systems strongly driven by stationary flow, we discuss an extended fluctuation- 
dissipation theorem (FDT). Beyond the linear response regime, the FDT for the stress acquires an 
additional contribution involving the observable that is conjugate to the strain rate with respect to 
the dissipation function. This extended FDT is evaluated both analytically for Rouse polymers and 
in numerical simulations for colloidal suspensions. More generally, our results suggest an extension 
of Onsager's regression principle to nonequilibrium steady states. 
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Introduction. - The fluctuation-dissipation theorem 
(FDT) holding for all systems slightly perturbed around 
their equilibrium state is one of the cornerstones of 
equilibrium statistical physics Specifically, for a 

system coupled to a heat bath at temperature T and 
setting Boltzmann's constant to unity throughout the 
manuscript, 

TR AJl (t -r) = (A(t)B(r)) = C A6 (t - r) (1) 

relates the response R A , h (£~ T ) = 3(A(t)) /5h(r) of an ob- 
servable A to a small perturbation h to equilibrium cor- 
relations. Crucially, these correlations involve the same 
observable A and the time-derivative B of another ob- 
servable that is conjugated to the perturbation h in the 
sense that upon a perturbation the energy U of the sys- 
tem transforms as U i— > U — Bh. In a stationary state, 
due to time-translational invariance both response and 
correlation function can only depend on the difference 
t — t > 0. The physical picture behind the FDT can be 
expressed by Onsager's regression principle: the decay 
of fluctuations created by a small external perturbation 
cannot be distinguished from the decay of spontaneous 
thermal fluctuations. 

Beyond the linear response regime, the FDT in the 
form ([T]) no longer holds. However, more than 30 years 
ago, Agarwal noted that any stationary Markov pro- 
cess obeys a generalized FDT [2] . This generalized FDT 
is obtained through linear response theory of stochastic 
processes where the probability distribution tp(t) obeys 
dtip = £ip with some time evolution operator C In 
the presence of a perturbation, one can decompose the 
evolution operator C — Cq + h5C into an unperturbed 
part Cq and a perturbation operator SC. The station- 
ary solution of the unperturbed system obeys Ctfij) s = 0. 
The generalized FDT then has the form of Eq. |T]) with 
B replaced by B* = Tip~ l 8C^) s . This form appears to 
be less useful than its equilibrium counterpart because 
the new conjugate variable B* has no independent phys- 
ically significant meaning. In order to equip the involved 
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variables with such a physical interpretation, we can take 
the effect of the driving into account as an additive ex- 
cess function. The general structure of such an extended 
FDT reads 

TR A , h (t) = (A(t)[B(0)-B(0)}) = C A z(t)-I A§ (t). (2) 

Interpreted in the spirit of Onsager's regression princi- 
ple, the FDT now states that the decay of forced fluc- 
tuations out of a nonequilibrium steady state cannot be 
distinguished from the decay of spontaneous fluctuations 
around B. Such a statement becomes physically signifi- 
cant only if the variable B carries a transparent physical 
meaning as an observable. One purpose of the present 
Rapid Communication is to show that for shear flow 
driven soft matter systems governed by stochastic dy- 
namics, the extended FDT |J2) acquires such a transpar- 
ent form: for a perturbation caused by a change of the 
strain rate the observable B becomes the stress that is 
conjugated to the strain rate with respect to the dissipa- 
tion function, the mean of which is related to the total 
entropy production. Such a characterization generalizes 
the identification of B as the local mean velocity in our 
previous study of a driven Langevin-type dynamics with 
h a small additional force [ij . A change of the frame of 
reference from the laboratory frame to the frame moving 
with this local mean velocity then restores the equilib- 
rium form of the FDT 0, i|. The extended FDT in its 
integrated form leads to an experimentally tested gener- 
alized Einstein relation Q. 

Our approach is complementary to the strategy of in- 
troducing an effective temperature to restore the equilib- 
rium form fT]) of the FDT even in nonequilibrium. The 
latter approach has been developed over the last decade 
for systems with a small heat flow into the reservoir corre- 
sponding to a small entropy production rate 0, H, 0, E3] • 
With such an effective temperature concepts from equi- 
librium statistical mechanics could be applied to driven 
systems even though a full microscopic understanding on 
the range of validity of this concept does not seem to 
have been reached yet. 

For the restricted but paradigmatic class of shear 
driven systems on which we will focus, quantitative 
progress has been achieved using the framework of mode- 
coupling theory. This includes the constitutive equa- 
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tion [Tl| using integration through transients [l2j and 
an FDT for the diffusion of a tagged particle [l3[ . Invari- 
ant quantities [3] constitute an exact result for systems 
driven through the boundaries with an unchanged bulk 
Hamiltonian, whereas in this Rapid Communication the 
system is driven through an imposed external flow. 

Soft matter under shear. - We consider soft matter 
systems such as colloidal suspensions or polymers that 
can be described as N interacting Brownian particles. 
The system is driven into a noncquilibrium steady state 
by shearing with a strain rate 7, resulting in an imposed 
flow profile u(r) = ^yG x of the solvent with unit vector 
e x . Picking out the ith particle, the force exerted by 
all other particles is Fj = —ViU, where the potential 
energy U(T) = J2(ij) u (\ r ij\) * s gi ven by the sum over all 
pairs (ij) with = — r--. The particles interact via 
an isotropic pair potential u(r) and the set of particle 
positions is denoted as T = {ri, . . . ,r^} [l9j . 

The response of an observable A(T) to a small, time- 
dependent variation of the strain rate is RA,-y{t — t;7), 
where the dependence on 7 emphasizes that such a re- 
sponse can be defined for any steady state, not only 
for 7 = corresponding to equilibrium. The mean 
(A(t)) = J dr A(T)ip{T,t) involves the time-dependent 
probability distribution ip(T,t). The perturbation oper- 
ator with respect to a small change of the strain rate is 
5C = — ^j^yid/dxi. In the linear response regime, the 
FDT 

TR AiJ {t-T;0) = (A(t)a xy (r)) (3) 

relates the response to a correlation function involving 
the stress 

A _^_ E aa2 ! M__£ >b . P) (4) 

m 1 131 i=i 

due to particle interactions. The resulting generalized 
fluctuation-dissipation relation reads 

TR An (t-r; 7 ) = (A(t)a* xy (T)) (5) 

with conjugate stochastic variable 

N 

a* y = T^-HCi> s = y t e x ■ TV, ln^ s . (6) 

i-l 

Applying the generalized Onsager principle by following 
Eq. ([2]), we split a* y — a xy — a xy into the stress (|4|) and 

N 

G X y = - Vt^x ■ Ft (7) 
i=l 

which involves the thermodynamic force T% = — Vj[Z7 + 
Tint's]- Both stresses have the same mean (cr xy ) = 
(cf xy ). In equilibrium, a xy = vanishes and hence Eq. (JSJ) 
reduces to Eq. §5§ as expected. 



We now provide for a xy a clear physical meaning con- 
necting it to the entropy production caused by the exter- 
nal flow u(r) . For overdamped dynamics, the dissipation 
function [151 ] 

1 N 

is related to the mean total entropy production rate 
through T(stot) = 2(W) 16]. In a nonequilibrium 
steady state, the sum (W) + A of mean dissipation func- 
tion and time-derivative of the "dynamical free energy" 
A = J dr ip[U + Tln-0] attains a minimum with re- 
spect to the local mean velocities {v^}, which then obey 
Vi = u(rj) + fio^i. A variation of W with respect to the 
strain rate leads to 

dW 1 A. du{r t ) 

-5- = V Vj - u(r;) • — = 0~ xy . 

07 mo ~[ 07 

In this sense, a xy is the variable conjugate to the strain 
rate 7 with respect to the dissipation function in analogy 
to B being the variable conjugate to h with respect to 
the energy. 

In the remainder of this Rapid Communication, we 
concentrate on response and correlation functions of A = 
B = (J xy - (o- xy ) and B = a xy — (a xy ), e.g., C(t) = 
(o~xy(t)a X y(0)) — (a xy ) 2 . To ease the notation we drop the 
subscripts denoting the observables. We first consider a 
Rouse polymer which allows for analytic expressions and 
then turn to numerical results for a colloidal suspension 
with nontrivial interactions. 

Rouse polymer. - Analytic expressions for correla- 
tion and response functions can be obtained for systems 
with quadratic interaction energies of the form U (T) = 
\ J2 a ^qIq w i tri amplitudes {q Q } of normal modes. The 
stress through interactions is o~ xy = ^ Q k a x a y a . Insert- 
ing the explicit form of ip B , Eq. j7]) becomes 

= H -1 2 (K 2 a x a y a + K a yl) . (9) 

A straightforward calculation [2FJ] based on the Smolu- 
chowski operator leads to a closed equation of motion for 
the correlation function 

C(t) = T 2 £[1 + «2 (3 + 2t/r a )]e- t ^ , (10) 

a 

where n a = jr a and r a = (2fi k a )^ 1 is half the re- 
laxation time of the corresponding mode. The response 
function R(t) = TJ2 a e~*/ r = is independent of the driv- 
ing. The excess I — C ~ TR oc j 2 is a quadratic function 
of the shear rate. While such quadratic behavior is ex- 
pected universally at small shear rates, for the Rouse 
polymer its persistence for large 7 depends on the Gaus- 
sian form of ip s . 

To obtain the universal expressions for a large num- 
ber of modes, we replace the summation by an inte- 
gration, ^2 a 1— > da, and set the relaxation times 
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FIG. 1: (color online) Integrated response \ vs. integrated 
correlation function K for a) Rouse polymer and b) a col- 
loidal suspension (see main text). Both functions have been 
normalized by Kaa = K(t — > oo). The straight solid line 
corresponds to the equilibrium FDT 



r a = Ti/a 2 , where the time scale is determined by the 
fundamental relaxation time t\. In addition, we consider 
integrated response \(t) = J* dr R(t) and correlation 
K{t) = J Q dr C(r). In Fig. 0Jl) the normalized inte- 
grated response x{K) is plotted against the value of the 
correlation function (parametrized by time) . If the equi- 
librium FDT © holds then this curve is a straight line. 
For increasing strain rates the deviation, and therefore 
the excess, increases. Moreover, for the Rouse polymer 
no regimes with constant slope corresponding to an ef- 
fective temperature can be discerned. 

Colloidal suspension. - We consider N colloidal par- 
ticles with diameter a suspended in a fluid. We assume 
the particles to interact through a repulsive, screened 
Coulomb pair potential 



u(r) = TZ 2 X B 



(l + na/2) 2 r ' 



(11) 



where Ab — 7 nm is the Bjerrum length in water at room 
temperature, Z is the effective surface charge, and k^ 1 
is the screening length. 

Response and correlation functions are obtained 
through simulations of a sheared dilute colloidal suspen- 
sion in a cubic box with side length L = 25a and vol- 
ume V = L 3 . We are interested in the bulk behavior 
and therefore we employ periodic Lees-Edwards bound- 
ary conditions in the simulation. The particle num- 
ber is N = 1000 corresponding to a volume fraction of 
<f> ~ 0.034. The screening length is set to k~ 1 = 0.15a and 
the effective surface charge is Z = 12000. The natural 
time scale r = a 2 /D is set by the time a particle needs 
to diffuse a distance equal to its diameter. To make con- 
tact with physical units, we choose a = 1 /im. In Fig. [2^), 
the mean stress is shown for different Peclet numbers 770. 
The straight line indicates the linear response behavior. 
The deviation of the mean stress from this line for large 
strain rates corresponds to shear thinning of the suspen- 
Mean values and correlation functions are ob- 
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tained from single runs with constant strain rate. The 
response function is determined as R(t) = dt(o~ xy (t)) /e 
after a jump 7 1— > 7 + e of the strain rate. Fig. ^Bp) shows 



FIG. 2: (color online) a) Mean stress (u xy ) of a colloidal 
suspension vs. the strain rate 7. The solid line indicates 
the linear response regime, b) Mean stress in the vicinity of 
"/to — 0.4. The solid line is a fit indicating the linear response 
of the nonequilibrium steady state. The error bars have been 
obtained as standard deviation by splitting the trajectory into 
eight segments, c) Pair distribution function g(r) for 7T0 = 1 
in the xy-p\&ne with z — 0. For parameters, see main text. 



the linear response of the nonequilibrium steady state 
corresponding to 7x0 = 0.4. In the simulation, we have 
chosen a step of £T$ = 0.01. 

The strong violation of the FDT can been seen in 
Fig. [3^i-c) , where response and correlation functions are 
plotted as functions of time for equilibrium (7 = 0) and 
two driven nonequilibrium steady states. For increas- 
ing strain rate, the deviation of the correlation function 
C(t) from the response function R(t) becomes larger, too. 
In Fig. [lb), the normalized integrated response x{K) IS 
shown as function of the integrated correlations. The 
small deviation between response and correlation func- 
tion in Fig. [3^,) is responsible for the fact that in Fig.QJ>) 
the equilibrium curve for 7 = lies slightly above the 
expected straight slope. Overall, the structure is compa- 
rable to the Rouse polymer in Fig.QJ,). Following Ref. [7}, 
one might even be tempted to identify a linear slope in 
the intermediate range corresponding to an effective tem- 
perature even though this is not the focus of the present 
work. 

Approximate excess function. - So far, we have ob- 
tained the excess as difference I = C — TR. In principle, 
we need the complete distribution ^s(r) to determine 
I independently using Eq. {J}. Such complete informa- 
tion is, however, neither available experimentally nor in 
computer simulations. Hence, approximate schemes will 
become important in future applications of the extended 
FDT. As a first step, we discuss here an approximation 
to a xy . Since the suspension is homogeneous, the one- 
point density pW(ri) = p = N/V is constant and the 
two-point density 

P (2) (n,r 2 ) =N(N-1) fdr 3 ---dr N </> 8 (r) = P 2 g(r) 



becomes a function of the displacement r = 17 — r 2 only. 
The factor N(N — 1) accounts for the possible permu- 
tations of the identical particles. The pair distribution 
g(r) as obtained from the simulation for the parameters 
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FIG. 3: (color online) Correlation function C(t), response function R(t), and approximated response function R(t) over time for 
a) equilibrium and b-c) two different nonequilibrium steady states, d) Comparison of integrated excess (through I — C — TR) 
and integrated approximated excess [using Eq. (|12[) ] vs. strain rate. 



introduced above is shown in Fig. [2];) (for hard spheres, 
cf. Ref. [3). 

We approximate the stationary distribution as ip s (F) = 
cxp{— u^2,uj\ w(rij)/T} with potential of mean force 
w(r) = — rin<7(r) [17(. This approximation effectively 
factorizes the probability distribution by using the cor- 
rect pair correlations and neglecting correlations between 
three and more particles, ft is motivated by the fact that 
the stress is determined by pair interactions only. In the 
parameter range we have studied we found a ~ 0.5. In- 
serting the pair approximation into Eq. ([6]) leads to 

The constant offset is adjusted for every strain rate such 
that (a* y ) = 0. We also employ a cut-off taking into 
account only neighboring particles within the first shell. 
Using this approximation, we can then both calculate the 
response function through ([5]) and access the stress (J7|) in 
the simulation. In Fig. [3K-c) the approximated response 
function R(t) is shown together with R(t) and the corre- 
lation functions. In Fig. [3jl), both the integrated excess 
J °° dr J(r) as well as the integrated approximate excess 
based on Eq. (|12|l are shown for different strain rates. 
For moderate to large strain rates, the pair approxima- 
tion works quite well. In the limit of vanishing strain 
rate, even though the volume fraction of the colloidal 
particles is low, the potential of mean force still deviates 
from the pair potential resulting in a break-down of this 



type of approximation for equilibrium. 

Concluding perspective. - For strongly driven soft mat- 
ter systems, we have discussed an extended FDT. Beyond 
the analytical and numerical data for two case studies, 
our general insight is twofold. First, beyond the linear 
response regime, the FDT acquires an additive contribu- 
tion which involves the stress that in the dissipation func- 
tion is conjugate to the strain rate. This result suggests 
more generally that the nonequilibrium form of the FDT 
involves the observable that is conjugate to the perturba- 
tion in the dissipation function. Such a scheme could be 
analogous to the pairing of observables conjugate with 
respect to energy in the equilibrium form of the FDT. 
Second, the additive contribution allows an interpreta- 
tion in the spirit of Onsager's regression principle: The 
decay of a spontaneous fluctuation around a nonequilib- 
rium steady state cannot be distinguished from the decay 
of a fluctuation forced by a small external perturbation. 
Whether these observations can be generalized to an even 
larger class of nonequilibrium systems remains to be in- 
vestigated both by further case studies and, more ambi- 
tiously, by an attempt to formulate a more formal theory 
for nonequilibrium steady states along these lines. 
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